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Approximate stabilization of an infinite dimensional quantum stochastic 

system 

Ram Somaraju, Mazyar Mirrahimi and Pierre Rouchon 



Abstract — We propose a feedback scheme for preparation 
of photon number states in a microwave cavity. Quantum 
Non-Demolition (QND) measurements of the cavity field and 
a control signal consisting of a microwave pulse injected into 
the cavity are used to drive the system towards a desired target 
photon number state. Unlike previous work, we do not use the 
Galerkin approximation of truncating the infinite-dimensional 
system Hilbert space into a finite-dimensional subspace. We 
use an (unbounded) strict Lyapunov function and prove that 
a feedback scheme that minimizes the expectation value of the 
Lyapunov function at each time step stabilizes the system at the 
desired photon number state with (a pre-specified) arbitrarily 
high probability. Simulations of this scheme demonstrate that 
we improve the performance of the controller by reducing 
"leakage" to high photon numbers. 



I. Introduction 

Quantum Non-Demolition (QND) measurements have 
been used to detect and/or produce highly non-classical states 
of light in trapped super-conducting cavities [1], [2], [3] 
(see [4, Ch. 5] for a description of such quantum electro- 
dynamical systems and [5] for detailed physical models with 
QND measures of light using atoms). In this paper we ex- 
amine the feedback stabilization of such experimental setups 
near a pre-specified target photon number state. Such photon 
number states, with a precisely defined number of photons, 
are highly non-classical and have potential applications in 
quantum information and computation. 

The state of the cavity may be described on a Fock space 
T-L, which is a particular type of Hilbert space that is used 
to describe the dynamics of a quantum harmonic oscillator 
(see e.g. [4, Sec 3.1]). The cannonical orthonormal basis 
for this Hilbert space consists of the set of Fock states 
{|0} , |1) , 1 2) , . . .}. Physically, the state |n) corresponds to a 
cavity state with precisely n photons. In this paper we study 
the possibility of driving the state of the system to some pre- 
specified target state |n). The feedback scheme uses the so 
called measurement back action and a control signal, which 
is a coherent light pulse injected into the cavity, to stabilize 
the system at the target state with high probability. 

Ram Somaraju and Mazyar Mirrahimi are with INRIA Rocquencourt, 
Domaine de Voluceau, B.R 105, 78153 Le Chesnay cedex, France, 
(ram. somaraju, mazyar.mirrahimi)@inria.fr 

Ram Somaraju and Mazyar Mirrahimi acknowledge support from 
"Agence Nationale de la Recherche" (ANR), Projet Jeunes Chercheurs 
EPOQ2 number ANR-09-JCJC-0070. 

P. Rouchon is with Mines ParisTech, Centre Automatique et Systemes, 
Mathematiques et Systemes, 60 Bd Saint Michel, 75272 Paris cedex 06, 
France, pierre.rouchon@mines-paristech.fr 

Pierre Rouchon acknowledges support from ANR (CQUID). 



Such feedback schemes for this experimental setup were 
examined previously in [6], [7]. The overall control structure 
used in [6] is a quantum adaptation of the observer/controller 
structure widely used for classical systems (see, e.g. [8, Ch. 
4]). The observer part consists of a discrete-time quantum 
filter, based on the observed detector clicks, to estimate the 
quantum-state of the cavity field. This estimated state is 
then used in a state-feedback based on Lyapunov design, 
the controller part. 

As the Hilbert space H is infinite dimensional it is difficult 
to design feedback controllers to drive the system towards 
a target state (because closed and bounded subsets of W 
are not compact). In [6], the controller was designed by 
approximating the underlying Hilbert space T~L with a finite- 
dimensional Galerkin approximation HN ma:c - Here, TiN max 
is the linear subspace of T~L spanned by the basis vectors 
|0) , |1) , .. . , \N max ) and N max > n, our target sate. Phys- 
ically this assumption leads to an artificial bound N max on 
the maximum number of photons that may be inside the 
cavity. In this paper we wish to design a controller for the 
full Hilbert space T-L without using the finite dimensional 
approximation. The need to consider the full Hilbert space 
is motivated by simulations (see Section |IV| > which indicate 
that using the controller designed on a finite dimensional 
approximation results in "leakage" to higher photon numbers 
with some finite probability. 

Controlling infinite dimensional quantum systems have 
previously been examined in the deterministic setting with- 
out measurements. Various approaches have been used to 
overcome the non-compactness of closed and bounded sets. 
One approach consists of proving approximate convergence 
results which show convergence to a neighborhood of the 
target state [9], [10]. Alternatively, one examines weak 
convergence for example, in [11]. Other approaches such as 
using strict Lyapunov functions or strong convergence under 
restrictions on possible trajectories to compact sets have also 
been used in the context of infinite dimensional state-space 
for example in [12], [13]. 

The situation in our paper is different in the sense that the 
system under consideration is inherently stochastic due to 
quantum measurements. The system may be described using 
a discrete time Markov process on the set of unit vectors 
in the system Hilbert space as explained in Section [IT] We 
use a strict Lyapunov function that restricts the system tra- 
jectories with high probability to compact sets as explained 



in Section III We use the properties of weak-convergence 
of measures to show approximate convergence (i.e. with 
probability of convergence approaching one) of the discrete 
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Fig. 1. The microwave cavity QED setup with its feedback scheme (in 
green). 



time Markov process towards the target state. 

We use a similar overall feedback scheme that is used 
in [6]. The entire feedback system is split into an observer 
part, a quantum filter, and a controller part based on a 
Lyapunov function. The quantum filter used to estimate the 
state is identical to the one used in [6] and we do not discuss 
the filter further in this paper. However we do not use the 
Galerkin approximation to design the controller. We show 
in Theorem 3.2 that given any e > 0, we can drive our 



system to the target state n with probability greater than 
1— e. Simulations (see Section[lV} indicate that this controller 
provides improved performance with lower probability of 
having trajectory escaping towards infinite photon numbers. 
The precise choice of Lyapunov function is motivated by [14] 
that uses a similar form of the Lyapunov function in a finite 
dimensional setting. 

A. Outline 

The remainder of the paper is organised as follows: in the 
following Section we describe the experimental setup and the 
Markovian jump dynamics of the system state. In Section III 
we state the main result of our paper including an outline of 



the proof of Theorem 3.2 We then present our simulation 
results in Section [IV] and then our conclusions in the final 
Section. 

II. System description 

The system, illustrated in Figure [T] consists of 1) a high- 
Q microwave cavity C, 2) an atom source B that produces 
Rydberg atoms, 3) two low-Q Ramsey cavities Ri and i? 2 , 
4) an atom detector D and 5) a microwave source S. The 
system may be modeled by a discrete-time Markov process, 
which takes into account the backaction of the measurement 
process (see e.g. [4, Ch. 4] and [6]). 

Rydberg Atoms are sent from B, interact with the cavity 
C, entangling the state of the atom with that of the cavity 
and are then detected in D. Each time-step, indexed by the 
integer k, corresponds to atom number k crossing the cavity 
and interacting with the cavity. The state of the cavity in 



time step k is described by a unit vector \ipk) € B\ for 
k = 1, 2, . . .. Here, B x = €H:\\ || = 1} is the set 
of possible cavity states. The change of the cavity state \ipk) 
at time-step k to the state |V>fe+i) at time-step k + 1 consists 
of two parts corresponding to the projective measurement of 
the cavity state, by detecting the state of the Rydberg atom 
in detector D and also due to an appropriate coherent pulse 
(the control) injected into C. 

Let a and a) be the photon annihilation and creation 
operators where a \n) — yfn \n — 1) and a) is the Hermition 
conjugate of a. Also, let N = a) a be the diagonal number 
operator satisfying N \n) = n\n). Let D a = exp(a(a^ — a)) 
be the displacement operator which is a unitary operator 
that corresponds to the input of a coherent control field of 
amplitude a that is injected into the cavity. The amplitude a 
of the coherent field is the control that is used to manipulate 
the system. Let M g = cos(0+ N<j>) and M e = sm(6+N(j)) be 
the measurement operators, where 9 and 4> are experimental 
parameters. Physically, the measurement operator M s , s € 
{e,g} correspond to the state of the detected atom in either 
the ground state \g) or the excited state |e) . 

We model these dynamics by a Markov process 



IV'fe+1/2) 
|^fc+i) 



\\M S |Vfe) II 
= D ak IVfc+1/2) • 



withprob. ||A/ S |^ fc )l| 2 (1) 
(2) 



Here s £ {e, g} and the control € R. 

Remark 2.1: The time evolution from the step k to k + 1, 
consists of two types of evolutions: a projective measurement 
by the operators M s and a coherent injection involving 
operator D a . For the sake of simplicity, we will use the 
notation of IV'fe+1/2) t0 illustrate this intermediate step. 

Remark 2.2: Let ftA\ be the set of all probability mea- 
sures on B\. Then the Equations <(TJ and |2]) determine 
a stochastic flow in A4% and we denote by Tk(no) the 
probability distribution of \ipk), given no, the probability 
distribution of |-0o)- 

III. Global (approximate) feedback stabilization 

We wish to use the control to drive the system into a 
pre-specified target state \n) with high probability. That is, 
we wish to show that the sequence converges to the 

set of probability measures flrx, where for all /i^ € fl^, 
Mood"-)) is big. 

In order to achieve this we use a Lyapunov function |5]) 
and at each time step k we choose the feedback control 
afc to minimize the Lyapunov function. Before discussing 
the choice of the Lyapunov function in Subsection |III-B| we 
recall some facts concerning the convergence of probability 
measures 



A. Convergence of probability measures 

We refer the interested reader to [15], [16] for results per- 
taining to convergence of probability measures. We denote 
by C the set of all continuous bounded functions on B\. 
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Definition 3.1: We say that a sequence of probability mea- 
sure {//„}^Lj C Mi converges (weak-*) to a probability 
measure /x G Mi if for all / G C 

lim E M „[/] =E M [/] 

n— >oo 

and we write 

It can be shown that if //„ ^-s- /ioo then for all open sets VT, 



limmi fx n (W) >Moo(W). 



(3) 



A set of probability measures S* C Mi is said to be tight [16, 
p. 9] if for all e > there exists a compact set K e C Bi 
such that for all fi E S, n(K e ) > 1 — e. 

Theorem 3.1 (Prohorov's theorem): Any tight sequence 
of probability measures has a (weak-*) converging subse- 
quence. 

We also recall Doob's inequality. Let X n be a Markov 
process on some state space X. Suppose that there is a 
non-negative function V{x) satisfying E[y(X 1 )|X = x)] — 
V(x) < 0, then Doob's inequality states 



supV(X n )> 1 \X = x) < 

n>0 



V(x) 
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(4) 



B. Lyapunov function and control signal 

We now introduce our Lyapunov function V and explain 
the intuition behind this peculiar form of this function. The 
function, V : Bi — > [0, oo] is defined 

oo 

V{\i>)) = ^a n K^|n>| 2 +5(co S 4 (^) + sin 4 (^)) 
-5(\\M g m\ 4 + \\M e \i:)\\ 4 ). (5) 



Here 



n<p, 



5 > is a small positive number and 



1 y,fi 
8 + 2^k= 


1 1 
-lk k 2 > 


if 7i = 




1 1 

k k 2 ' 


if 1 < < n 


o, 




if n = n 




fe ^ fe 2 ' 


if n > n 



(6) 



We set Z>(V) C Bi to be the set of all \tp) £ B x where the 
above Lyapunov function is finite. We note that coherent 
states, which are states that are of relevance in practical 
experiments are in T>(V). 

We choose a feedback that minimizes the expectation 
value of the Lyapunov function in every time-step k. Indeed, 
applying the result of the fc'th measurement, we know the 
state Vfc+1/2) and we choose a k as follows 



a k = argmin V (D a ^fc+i/a}) 



(7) 



for some positive constant a. 

Remark 3.1: The Lyapunov function is chosen to be this 
specific form to serve three purposes - 

1) We choose the sequence er„ — > 00 as n — > 00. 
This guarantees that if we choose a k to minimize the 



expectation value of the Lyapunov function then the 
trajectories of the Markov process are restricted to a 
compact set in Bi with probability arbitrarily close to 
1. This implies that the w-limit set of the process is 



non-empty (see Step 2 in the Proof of Theorem 3.2 1. 
2) The term -S(\\M g \ip) || 4 + \\M g \tp) || 4 ) is chosen 
such that the Lyapunov function is a strict Lyapunov 
functions for the Fock states. This implies that the 
support of the w-limit set only contains Fock states 



(see Step 3 in the Proof of Theorem 3.2 1. 
3) The relative magnitudes of the coefficients u n have 
been chosen such that V^(|n)) is a strict global mini- 
mum of V. Moreover given any M > n we can choose 
5, a such that for all M > m ^ n, and for all in 
a neighborhood of |m), V(D a does not have a 
local minimum at a = 0. This implies that if \ip k ) 
is in this neighborhood of |m) then we can choose 
an at € [—a, a] to decrease the Lyapunov function 
and move \ipk) away from \m) by some finite distance 
with probability 1 (see Steps 4 and 5 in the Proof of 
Theorem 1 3. 2b. 



C. Main Result 

We make the following assumption. 

Al The eigenvalues of M g and M e are non-degenerate. This 
is equivalent to the assumption that 7r/0 is not a rational 
number. 

The quantum filter uses the statistics of the measurement of 
whether the atom is in the ground or excited state to estimate 
the cavity's state. Therefore if one of the eigenvalues of M g 
(or M e ) is degenerate then the measurement statistics will be 
the same for more than one photon number state. Therefore 
it is not possible to control the system effectively in this 



case (However, as explained in Remark 3.2 below, we may 
weaken this assumption slightly). 

The following Theorem is our main result. 

Theorem 3.2: If we assume Al to be true then given any 
e > and C > 0, there exist constants 6 > and a such 
that for all /j satisfying E M [V] < C, T n (fi) converges to a 
limit set Q. Moreover for all fi^ € f2, \ip) £ supp(/Xoo) only 
if \ijj) is one of the Fock states \n) and 

Mco({|rl>}) >l-e. 
The proof is split into 5 steps: 

1) V(\ipk)) is a super-martingale that is bounded from 
below. 

2) The sequence of measures T k (/x) is tight and therefore 
has a converging subsequence. Hence the set il is non- 
empty. 

3) If Tk l (/i) — > /ioo tnen the support set of fioo only 
consists of Fock states. 

4) Let M', C > be given. Then for all W > m ^ n, 
5 and a may be chosen small enough such that for 
k > small enough and all in the neighborhood 

K - {IV^> : II IV)-H II < > v(H)-k} 

(8) 
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of |m), satisfying V(\ip)) < C, we have for \a\ < a 
the polynomial approximation 



V(D a |^)) = + 0(a 3 ) + 0(5) 

and f 2 (\ijj)) < 7 < for some constant 7. The term 
0(a 3 ) only depends on C and not on \ip) and the 
term 0(5) is independent of both and C". 
5) Because 7 is negative, we can choose a and (5 small 
enough such that the probability of convergence to the 
Fock states |m) for m 7^ n may be made arbitrarily 
small. Therefore 



Mood")) 



1 



Mood™)) 



m=0 



may be made arbitrarily big. 
Below we sketch the proofs of each of the above steps. 
The interested reader is referred to [17] for further details 
on the proof which are beyond the scope of a short note. 
Proof: [Proof of step 1] We can write 

where, 



min E [V (D a (|^ k+ i/2») \ \tpk}] 

a£[-a,Q] 



-E[V(|D (V fc+ i /a »||V*)] 
tfaflVfc)) = E[y(D (|V* + i/a>)) | ItM] 



(9) 



It is obvious that Ki(\ip)) < and after simple but tedious 
manipulations, we get 

im\2 



-2(||M 9 2 |^)|p 



Tr{M g V}Tr{M|p} 



< 0. (10) 



Therefore, y(^fe) is a super-martingale. ■ 
Proof: [Proof of step 2] Let e > be given. Because 
V(\ipk)) is a supermartingale, Doob's inequality (JTTJ gives 
us 

8UpV(|^ fc ))>^ 



< e. 



(ID 



If we set, 



K e = {\iP):V(m<E ll [V}/e}) 



then for all k > 0, [rfc(/i)](iQ > 1 — e. Because, the 
sequence er n — > 00 as n — >• 00, the set X e can be shown 
to be pre-compact in W. We can now apply Prohorov's The- 
orem 3.1 to show that r n (/j,) has a converging subsequence. 
Therefore the limit set ft = {/loo e Mi : T^^t) /Woo} is 
non-empty. ■ 
Proof: [Proof of step 3] Suppose some subsequence 
of Tfc(/i) converges to /ioo £ f2. From step 1 we have 
^'i(IV'fc)) + K 2 (\ipk)) — > as fc — > 00 and because K\ 
and K 2 are both non-negative we have 

lim E rkW [K 2 ] = 0. 



But, from ( |T0] > and the boundedness of M g and M e , we 
know that K 2 is a continuous function on "H. Therefore from 



Definition 3.1 of (weak-*) convergence of measures we get 



E, 



[K 2 ] = 0. 



But K 2 {\i>)) = implies \\M 2 g \^)\\ 2 = \\M g ft,) 
Cauchy- Schwartz inequality gives 



(12) 
I 4 . The 



\M 2 q |V) 



\M 2 g \^)\\ 2 \ 



> 



l^)H 2 

<v>M g 2 |M g v) • m) 
\im^)\ 2 

l|M 9 |^)|| 4 . 



with equality if and only if and M g \ip) are co-linear. 
Therefore K 2 (\i/j)) = implies (by Assumption^!.) that 
is a Fock state. Hence from ( fT~2] > we can conclude that the 
support set of fi^ only consists of the set of Fock states. ■ 
Proof: [Proof of step 4] 
Set 

OO 

V(m^J2^\(D a ^\m)\ 2 

n=0 

It can be shown [17] that V(D a is an analytic function 
of a if satisfies V(\ip)) < 00. Moreover, for all 
satisfying V(\ip)) < C we have the second order polynomial 
approximation 



V(D a |^)) = 



2 i 

Ea 
If 



t V(D a \i;))\ a=0 + O(a 3 ) 



for all \a\ < a. In particular the O(a) term only depends 
on C and is independent of Here V^(-)| Q= o is the i th 
derivative of (•) w.r.t. a evaluated at a = 0. 

If we let \ip) = Y^=Q c n\ n ) an d recall that D a = 
exp(a(a — a')) then after some manipulations, we get 



v 2 v(D a m\ 



a=0 
2 / 



= 2^ \ c n\ 2 ((n + 1)^+1 + na n -i - (2n + l)a n ) 
n=0 

+ Rc{c n -iCn +1 }y/n(n + 1)(<t„-i + (J n +i - 2ct„). 
If n 7^ n and n > 2 we have 

(n + l)cr„+i + ucr n _i - (2n + l)a r 

and for n = 0, 1 we get 



-1 
n(n + 1) 



(n + l)(T n+ i + na n -i - (2n + 1)ct„ = 



-1 



For any Fock state \m) with m 7^ n, c n = S mn , 
S mn is the Kronecker-delta function and we have 

7 2?~- • - 1 



where 



V^(i) Q |m))| Q=0 = - 



< 0. 



m(m + 1) 

Because the terms J2 n l c ™l 2 and En Re { c "+i c *i-i} are 
bounded by the || • ||-norm in T~L, it can be shown that for 



k small enough we have V a V{D a |V , ))| Q ,- < 
in the neighborhood V? n of \m), where Vf n is given as in 
Equation ([8]). 
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But, 

V%V(D a m\ a=Q = V 2 a V(D a m\ a=Q + 0(5). 

Hence, given any M > h, step 4 above is true with 7 = 
1 ■ 

2M(M+1)- 

Proof: [Proof of step 5] Let e > be given. We 
show that ^oo({|")}) > 1 — e. From step 3 we know that 
the support of fi^ only consists of Fock states. Therefore 
using ([3]), we only need to show that there exists an open 
neighborhood W of {\m) : m ^ n} such that for k big 
enough the [T k (p)]{W) < e. 

We construct the set W using two disjoint parts W\ 
and W%. We first show that there exists a M big enough 
and a neighborhood W\ of {\M) , \M + 1) , . . .} such that 
[rfe( y u)](W / i) < e/2 for all k. We then construct a neigh- 
borhood W2 of {|m) : < m < M,m 7^ n} such that 
[r fc (/i)](W 2 ) < e/2 for /c large enough. 

a) Construction of W%: Because a m — > 00 there 
exists an Af large enough such that for all to > M, 
cr m > We can choose a small enough neighborhood 
W x of \\M) , \M + 1) , . . .} such that for all |^) in this 
neighborhood, 

\wn - 2 - e ^ 2 

Because E^IV] < C, Doob's inequality implies the proba- 
bility of V(\ip k )) > C/(e/2) is less than e/2. Therefore, 



[r fc (/i)](^i) < 



(13) 



Construction of W%: We show that for k small 
enough we can choose 



M-l 



w 2 = (J K 



m=0 

where is as in ([8]). 

From Doob's inequality, we have 

F*Ml({W:W))>4})<£/2. (14) 

for all fc. Therefore we can complete the proof if we show 
that for n small enough 



lim [r fe ( M )](v«)=o, 



where 



v K = v K n 



e72 



In Step 4 we set C" = C/(e/2) and M' = M and let 
k be small enough so that V? n is as given in step 4. Then, 
because 7 < 0, we can choose a and 8 small enough so that 
there exists a constant c > such that for all £ V^, 
V(D a \tp}) — V(\ip)) < — c, for some a € [—a, a]. Because 
M g and M e are bounded operators and M g \m) = |m) and 
M e |m) = I to), k' can be chosen small enough such that if 
\ii k ) € V<i then l^fe+1/2) G with probability 1. 



We claim that r /c (/i)(V i 
because for all € 
1 and this implies 



— > 0. To see this note that 
1/2) G Vm w i m probability 



FdV'fc+i}) - V^d^fc)) < -c with probability 1. 

So if \tjjk) G Vm' 2 then the Markov process is outside the 
set within a finite number of steps less than [C"/c] 
with probability 1. So if /ifc(Vm' 2 ) does not approach zero, 

A /2 

then the Markov process must enter the set V m from 
outside the set V m infinitely many times. But by Doob's 
inequality (JTTJ the probability of this happening once is less 
than (V(|m)) - n')/{V{\m)) - n'/2) < 1. Therefore the 
probability of this happening infinitely many times is zero. 
Thus lim / t- i .oo[rA : (^)](Vm' 2 ) = 0. This combined with ( p"3j ) 
and ( fT4| i gives, /.ioo{\m) ■ m ^ n} < e. 
Therefore, 



|n)) = l- ^Mcx>(|m))>l-e. 

m— 



Remark 3.2: In step 2 we show that the only vectors in 
the support of /i^ are those corresponding to eigenvector 
of M s . We then used assumption Al to claim that the only 
eigenvectors of M s are the Fock states. We can however 
weaken this assumption to the following: eigenvalues corre- 
sponding to eigenvectors |m), to < M are non-degenerate. 
This is because, we can show that if some eigenvector 
is in the span of the set {\M) , \M + 1) , . . .} then using the 
same argument as that used for |m) , m > M, we can show 
that the probability of is small. This is significant for 
cases where M g is a more complicated non-linear function 
of AT, as is the case in a practical system. 

IV. Simulations 



To illustrate Theorem 3.2 we performed closed-loop simu- 



lations of the controller designed using the finite-dimensional 
approximation [6] and the one in Theorem |3.2| Both simu- 
lations were performed on a system truncated to 21 photons. 
However the quantum filter (and therefore the controller) was 
truncated to 10 photons. 

The initial state was chosen to be the coherent state having 
an average of n = 3 photons: 



I'M 



The measurement 

cos (y/2(N -n)/5 + 



We take 



operators 

M P = sin I 



are M g = 
(V2(JV-n)/5+f). 
a = ± and 6 = (1/10(10 + l))/2 to ensure 
the Lyapunov function is strictly concave near the Fock 
states I to), to 7^ n. To compute the feedback law given 
by the minimisation (7}, we approximate, for each step k, 
[—a,+a] 3at4E \ V ({tpk+i)) \ \~>pk}] by the polynomial 
of degree two with the same first and second order 
derivatives at a k = 0. Figure [2] shows good convergence 
properties of such feedback strategy with an average 
asymptotic value of | {tt.| -i/?) 1 2 close to 1. The remaining 
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Fidelity between |w > and the goal Fock state 




Fig. 2. Simulation with a truncation to 20 photons of the system and 9 
photons of the filter for the feedback law jTJ; in blue|(?i|i/>£,}| 2 (n = 3) for 
each realization ; in red average over the 100 realizations of |(»2|i/>j.}| 2 . 

Fidelity between \y k ) and the goal Fock state 




200 250 300 

Step number k 



trajectories that do not converge to \n) can be interpreted as 
the e in theorem [3. 21 

Figure [3] is devoted to similar simulations but with the 
feedback law of [6], [7] based on a finite dimensional model: 



< A: 



Ctk 



a if |(n|-0fc+i/2>| 10 , 

[i>k+i/2 | [l"> <n| ,a i -a,]ip k+1/2 ^ 
4n+2 



otherwise. 



(15) 



The average asymptotic value of | (n\ip) | is then around 0.95 
with this "finite dimensional" feedback. Around 5% of the 
trajectories do not converge towards |n) and escape towards 
high photon numbers. Figure |4] shows a typical example of 
such a trajectory which converges towards photon number 
15 and 20. 

V. Conclusion 

In this paper we examine the stabilization of a quantum 
optical cavity at a pre-specified photon number state |n). In 
contrast with previous work, we designed a Lyapunov func- 
tion on the entire infinite dimensional Hilbert space instead 
of using a truncation approximation. The Lyapunov function 
was chosen so that it is a strict Lyapunov function for the 
target state and the feedback consisted of a control that 
minimizes the expectation value of the Lyapunov function 
at each time-step. Simulations indicate that this feedback 
controller performs better than the one designed using the 
finite dimensional approximation. 



Fig. 3. Simulation with a truncation to 20 photons of the system and 
9 photons of the filter for the "finite dimensional" feedback law (15) ; in 
blue|(n|i/>fc}| 2 (n = 3) for each realization ; in red average over the 100 
realizations of |(n|?/>t,)| 2 . 
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Fig. 4. An example of a trajectory of the finite-dimensional controller 
demonstrating escape to high photon numbers. 
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